Relation between baryon number fluctuations and experimentaUy observed proton 
number fluctuations in relativistic heavy ion coflisions 
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We explore the relation between proton and nucleon number fluctuations in the final state in rela- 
tivistic heavy ion collisions. It is shown that the correlations between the isospins of nucleons in the 
final state are almost negligible over a wide range of collision energy. This leads to a factorization 
of the distribution function of the proton, neutron, and their antiparticles in the final state with 
binomial distribution functions. Using the factorization, we derive formulas to determine nucleon 
number cumulants, which are not direct experimental observables, from proton number fluctua- 
tions, which are experimentally observable in event-by-event analyses. With a simple treatment for 
strange baryons, the nucleon number cumulants are further promoted to the baryon number ones. 
Experimental determination of the baryon number cumulants makes it possible to compare various 
theoretical studies on them directly with experiments. Effects of nonzero isospin density on this 
formula are addressed quantitatively. It is shown that the effects are well suppressed over a wide 
energy range. 
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I. INTRODUCTION 

Now that the observation of the quark-gluon matter 
in relativistic heavy ion collisions is established for small 
baryon chemical potential (/is) ^ challenging exper- 
imental subject following this achievement is to reveal 
the global structure of the QCD phase diagram on the 
temperature (T) and /zb plane. In particular, finding the 
QCD critical point(s), whose existence is predicted by 
various theoretical studies 0, Q, is one of the most in- 
triguing problems. Since the hb of the hot medium cre- 
ated by heavy-ion collisions can be controlled by varying 
the collision energy per nucleon pair, ^snn, the /ie de- 
pendence of the nature of QCD phase transition should 
be observed as the ^/suN dependence of observables. An 
experimental project to explore such signals in the en- 
ergy range lOGeV < ^/suN < 200GeV, which is called 
the energy scan program, is now ongoing at the Rela- 
tivistic Heavy Ion Collider (RHIC) [JQ. Experimental 
data which will be obtained in future experimental facil- 
ities designed for lower beam-energy collisions will also 
provide important information on this subject @. 

Observables which are suitable to analyze bulk proper- 
ties of the matter around the phase boundary of QCD in 
heavy ion collisions are fluctuations Q- Experimentally, 
fluctuations are measured through event-by-event analy- 
ses j^]. Theoretically, it is predicted that some of them, 
including higher-order cumulants, are sensitive to criti- 
cal behavior near the QCD critical point [sl-fl^. and/or 
locations on the phase diagram, especially on which side 
the system is, the hadronic side or the quark-gluon side 

Mm. 



Among the fluctuation observables, those of conserved 
charges are believed to possess desirable properties to 
probe the phase structure in relativistic heavy ion col- 
lisions. One of the advantages of using the conserved 
charges is that the characteristic times for the variation 
of their local densities are longer than those for non- 
conserved ones, because the variation of the local den- 
sities of conserved charges are achieved only through dif- 
fusion ■ The fluctuations of the former thus can 
better reflect fluctuations generated in earlier stages of 
flreballs, when the rapidity coverage is taken sufficiently 
large. From a theoretical point of view, an important 
property of the conserved charges is that one can de- 
fine the operator of a conserved charge, Q, as a Noether 
current. Moreover, their higher-order cumulants, {SQ"')c, 
are directly related to the grand canonical partition func- 
tion Zifi) = Tre-''(^-'^'3) as 
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with H and fi being the hamiltonian and the chem- 
ical potential associated with Q, respectively. These 
properties make the analysis of cumulants of conserved 
charges well defined and feasible in a given theoreti- 
cal framework. For example, they can be measured 
in lattice QCD Monte Carlo simulations [20l - [2i |. The 
relation Eq. ([l} also provides an intuitive interpre- 
tation for the behavior of higher-order cumulants of 
conserved charges. For instance, the third-order cu- 
mulant of the net baryon number, Nj^'^^\ satisfies 
{iSN^^''^f), = r9(((5iV^""*))2),/a^B. This formula 

means that {{SN^'^*'^)^)^ changes its sign around the 
phase boundary on the T-^b plane where the baryon 
number susceptibility ((<5-/Vb"°*-*)^) has a peak structure 
[TtI I . The change of the sign of observables like this will 
serve as a clear experimental signal (T7l-[l9|. 
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QCD has several conserved charges, such as baryon 
and electric charge numbers and energy. Among these 
conserved charges, theoretical studies suggest that the 
cumulants of the baryon number have the most sensi- 
tive dependences on the phase transitions and phases of 
QCD. In order to see this feature, let us compare the 
baryon number cumulants with the electric charge ones. 
First, the baryon number fluctuations show the critical 
fluctuations associated with the QCD critical point more 
clearly. Although the baryon and electric charge num- 
ber fluctuations diverge with the same critical exponent 
around the critical point, it should be remembered that 
this does not mean similar clarity of signals for the crit- 
ical enhancement in experimental studies. Fluctuations 
near the critical point are generally separated into singu- 
lar and regular parts, and only the former diverges with 
the critical exponent. The singular part of the electric 
charge fluctuations is relatively suppressed compared to 
the baryon number ones, because the formers contain 
the isospin number fluctuations which are regular near 
the critical point Q . The additional regular contribution 
makes the experimental confirmation of the enhancement 
of the singular part difficult, and this tendency is more 
pronounced in higher-order cumulants [l7j . While it is 
known that the proton number fluctuations in the final 
state also reflect the critical enhancement near the crit- 
ical point ^9], as we will show later the baryon number 
fluctuations are superior to this observable, too, in the 
same sense. Second, the ratios of baryon number cu- 
mulants [l6j behave more sensitively to the difference 
of phases, i.e., hadrons, or quarks and gluons. This is 
because the ratios are dependent on the magnitude of 
charges carried by the quasi-particles composing the state 
[HQ [3, while the charge difference between hadrons 
and quarks is more prominent in the baryon number. 

Experimentally, however, the baryon number fluctu- 
ations are not directly observable, because chargeless 
baryons, such as neutrons, cannot be detected by most 
detectors. Proton number fluctuations can be measured 
[1, Q, and recently its cumulants have been compared 
with theoretical predictions for baryon number cumu- 
lants. Indeed, in the free hadron gas in equilibrium the 
baryon number cumulants are approximately twice the 
proton number ones, because the baryon number cumu- 
lants in free gas are simply given by the sum of those 
for all baryons, and the baryon number is dominated by 
proton and neutron numbers in the hadronic medium 
relevant to relativistic heavy ion collisions. In general, 
however, these cumulants behave differently. In fact, we 
will see later that the non-thermal effects which exists in 
baryon number cumulants are strongly suppressed in the 
proton number ones. 

In heavy ion collisions, because of the dynamical evo- 
lution the medium at kinetic freezeout is not completely 
in the thermal equilibrium. The original ideas to exploit 
fluctuation observables as probes of primordial properties 
of flreballs [l^, (TJ] are concerned with this non-thermal 
effect encoded in the flnal state as a hysteresis of the time 



evolution. To observe such effects, it is highly desirable 
to measure baryon number cumulants that is expected to 
retain more effects of the phase transition and the sin- 
gularity around the critical point. The experimental de- 
termination of baryon number cumulants also makes the 
comparison between experimental and theoretical stud- 
ies more robust, since many theoretical works including 
lattice QCD simulations are concerned with the baryon 
number cumulants, not the proton number ones. 

In Ref. [25)], the authors of the present paper have ar- 
gued that, whereas the baryon number cumulants are not 
the direct experimental observables as discussed above, 
they can be determined in experiments by only using the 
experimentally measured proton number fluctuations for 
-^SNN ^ lOGeV. The key idea is that isospins of nucleons 
in the flnal state are almost completely randomized and 
uncorrelated, because of reactions of nucleons with ther- 
mal pions in the hadronic stage, as will be elucidated 
in Sec. |lTl This leads to the conclusion that, when A^n 
nucleons exist in a phase space of the final state, the 
probability that Np nucleons among them are protons 
follows the binomial distribution. More generally, the 
probability distribution that Np protons, Nn neutrons, 
Np anti-protons, and Nn anti-neutrons are found in the 
final state in a phase space is factorized as 



rNiNp,Nn,Np,Nn) 
= TiN^,N^)BriNp; N^)BriNp; N^), 



(2) 



where the nuclcon and anti-nucleon numbers are An = 
Ap -I- A„ and A^jj = Np + Nn, respectively, and 

Br{k;n)^—^^r'^{l^rr-'' (3) 



k\in-k) 



is the binomial distribution function with probabilities 
r = {Np)/{Nys) and f = (A^p)/(A^n). The function 
J^{Nn, N-^) describes the distribution of nucleons and 
anti-nucleons and the correlation between them in the 
final state, which are determined by the dynamical his- 
tory of fireballs. Using the factorization Eq. ([2]), one can 
obtain formulas to represent the (anti-)nucleon number 
cumulants by the (anti-)proton number ones, and vice 
versa; whereas the neutron number is not determined 
by experiments, this missing information can be recon- 
structed with the knowledge for the distribution function, 
Eq. ([2]). The (anti-)nucleon number in Eq. ^ can fur- 
ther be promoted to the (anti-)baryon number in practi- 
cal analyses with a simple treatment for strange baryons 
to a good approximation. These formulas enable to de- 
termine the baryon number cumulants solely with the ex- 
perimentally measured proton number fiuctuations, and, 
as a result, to obtain insights into the present experimen- 
tal results on the proton number cumulants. 

The main purpose of the present paper is to elabo- 



rate the discussion in Ref. |25| with some extensions. 
In Ref. [1^ the formulas are derived only for isospin 
symmetric medium. In the present study we extend 
them to incorporate cases with nonzero isospin densi- 
ties. With the extended relations, it is shown that the 
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effect of nonzero isospin density is well suppressed for 
i7sNN ^ lOGeV. The procedures of the manipulations 



and discussions omitted in Ref. [25| are also addressed in 
detail. 

In the next Section, we show that the factorization 
Eq. ([2|) is well applied to the nucleon and baryon dis- 
tribution functions in the final state in heavy ion colli- 
sions. We then derive formulas to relate baryon and pro- 
ton number cumulants in Sec. IIIII In Sec.l IVl we discuss 
the recent experimental results at STAR [J, |5| using the 
results in Sec. IIII[ and possible extensions of our results. 
The final section is devoted to a short summary. 

Throughout this paper, we use Nx to represent the 
number of particles X leaving the system after each col- 
lision event, where X = p, n, N, and B represent pro- 
ton, neutron, nucleon, and baryon, respectively, and their 
anti-particles, p, n, N, and B. The net and total numbers 



are defined as N 
respectively. 



(not) 
X 



N- 



X 
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II. DISTRIBUTION FUNCTION FOR PROTON 
AND NEUTRON NUMBERS 

In this section, we discuss the time evolution of the 
proton and neutron number distributions in the hadronic 
medium generated by relativistic heavy ion collisions, 
and show that the nucleon distribution in the final state 
in a phase space is factorized as in Eq. ^ at sufficiently 
large ^/sKN■ In Sec. Ill A| as a preliminary example we 
show that Eq. ^ is applicable to the equilibrated free 
hadron gas in the ranges of T and /ib relevant to relativis- 
tic heavy ion collisions. We then extend the argument to 
the distribution function in the final state in relativistic 
heavy ion collisions in Sec. Ill Bl 



A. Free hadron gas in equilibrium 

Let us first consider nucleons in the free hadron gas in 
equilibrium. For T and /ie which are relevant to rela- 
tivistic heavy ion collisions, the nucleon mass mN satis- 
fies TTiN — l/iel 2> T. One thus can apply the Boltzmann 
approximation for the distribution functions of nucleons. 
The number of particles in a phase space, N, which obey 
Boltzmann statistics is given by the Poisson distribution. 



PxiN) = 



Nl 



(4) 



with the average A — (N) = X]jv -^^■^(-^)- Accordingly, 
the probability to find Np (Np) protons (anti-protons) 
and Nn (Nn) neutrons (anti- neutrons) in the phase space 
is given by the product of the Poisson distribution func- 
tions, 

VHG{Np,Nn,Np,Nn) 

= P^x,^ {Np)P^x„^ {N„)P^x,) (A^p)P(^«) (iV^). (5) 



The product of two Poisson distribution functions sat- 
isfies the identity, 

PxANi)PxAN2) 

= Px,+x,{Ni + N2)Bx,/ix,+x,){Ni;Ni + N2), (6) 

where Br{k\n) is the binomial distribution function 
Eq. (O. Using Eq. ©, Eq. (O is rewritten as 

VYiG{Np,Nn,Np,Nn) 

= P(N^){m)P(N^){N^)Br{Np-N^)B,{Np;N^), (7) 



where N] 



N 



N„ + Nr, and iV, 



N 



Np + Nn are the nu- 
cleon and anti-nucleon numbers, respectively, and r — 
{Np)/{Nt^) and f = {Np)/{N^). Equation ^ shows that 
the distribution of nucleons in the free hadron gas is fac- 
torized using binomial functions as in Eq. ([2]) with 



TiN^,N^) = P(^^,)(A^N)P(JVs)(^^N)■ 



(8) 



The appearance of the binomial distribution functions 
in Eq. ([7]) is understood as follows. When one finds a 
nucleon in the hadron gas, the probability that the nu- 
cleon is a proton is r. The isospins of all nucleons found 
in the phase space, moreover, are not correlated with one 
another as a consequence of Boltzmann statistics and the 
absence of interactions. Once A^n nucleons are found in 
the phase space, therefore, the probability that Np parti- 
cles are protons is a superposition of independent events 
with probability r, i.e., the binomial distribution. 

We note that the above discussion is not applicable 
when the condition ton — |mb| 3> T, required for Boltz- 
mann statistics, is not satisfied. When quantum corre- 
lations of nucleons arising from Fermi statistics are not 
negligible, the isospin of each nucleon can no longer be 
independent. As long as we are concerned with the range 
of T and fiB which can be realized by relativistic heavy 
ion collisions, however, the condition for the Boltzmann 
approximation is well satisfied except in very low energy 
collisions [26j |. 



B. Final state in heavy ion collisions 

Next, we consider the nucleon distribution functions in 
the final state in heavy ion collisions. We show that the 
nucleon distribution in this case is also factorized as in 
Eq. ([5]), by demonstrating that the isospins of all nucleons 
in the final state are random and uncorrelated. 



1. A(1232) resonance 

The key ingredient to obtain the factorization Eq. ([2]) 
in the final state in relativistic heavy ion collisions is Ntt 
reactions in the hadronic stage mediated by A(1232) res- 
onances having the isospin / = 3/2. As we will see later, 
this is the most dominant reaction of nucleons in the 
hadronic medium. This is because i) the cross section of 
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Ntt — )• A reactions exceeds 200nib = 20fni and is com- 
parable with NN and NN reactions for Piab ^ 300MeV 
|27| . and ii) the pion density dominates over those of all 
other particles in the ranges of T and fiB accessible with 
heavy ion collisions at ^snn ^ lOGeV; at the top RHIC 
energy, the density of pions is more than one order larger 
than that of nucleons. We shall show below that these 
reactions frequently take place even after chemical freeze- 
out in the hadronic medium during the time evolution of 
the fireballs. 

The Ntt reactions through A contain charge exchange 
reactions, which alter the isospin of the nucleon in the 
reaction. The reactions of a proton to form A are: 



P + TT ' 
P + TT 

p + n 







A++ -^p + n+, 

A+ ^p(n)+^°(7r+), 

A" ^ p{n)+TT-{TT°). 



(9) 
(10) 
(11) 



Among these reactions, Eqs. ((TU)) and pT|) are responsible 
for the change of the nucleon isospin. The ratio of the 
cross sections of a proton to form A++, A+, and A° is 3 : 
1:2, which is determined by the isospin SU(2) symmetry 
of the strong interaction. The isospin symmetry also tells 
us that the branching ratios of A+ (A") decaying into 
the final state having a proton and a neutron are 1 : 2 
(2 : 1). Using these ratios, one obtains the ratio of the 
probabilities that a proton in the hadron gas forms A+ 
or A° with a reaction with a thermal pion, and then 
decays into a proton and a neutron, respectively, Pp^p 



and Pri 



P, 



p^p 



P 



= 5:4, 



(12) 



provided that the hadronic medium is isospin symmetric 
and that the three isospin states of the pion are equally 
distributed in the medium. Because of the isospin sym- 
metry of the strong interaction one also obtains the same 
conclusion for neutron reactions: 



P, 



5 : 4. 



(13) 



Similar results are also obtained for ant i- nucleons. Equa- 
tions (IT2|) and (fT3|) show that these reactions act to ran- 
domize the isospin of nucleons during the hadronic stage. 

2. Mean time 

Next, let us estimate the mean time of these reactions. 
Assuming that pions are thermally distributed, the mean 
time TA of a nucleon at rest in the medium to undergo a 
reaction Eq. ([TU)) or PT|) is given by 



{2ny 



(14) 



with the Bose distribution function n{E) — {e-^^"^ ~1) ^, 
the pion momentum A:^, the pion velocity i;^ = kTr/Ej^, 
Ett = \/ m\ + fc^ , and the pion mass ttIt^. (T(i?c.m.) is the 
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FIG. 1: Mean time ta of a rest nucleon to form A+ or A" 
in the hadronic medium as a function of temperature T . 



sum of the cross sections for Ntt reactions producing A+ 
and A° for the center-of-mass energy iJc.m. = [(^N + 
i?^)2 — fcj]^/^ with the nucleon mass ttin- For the cross 
section CT(i?c.m.), we assume that the peak corresponding 
to A(1232) resonance is well reproduced by the Breit- 
Wigner form, 



f (£^c.m.) Ca 



rV4 



{E, 



mA) 



r2/4' 



(15) 



which is a sufficient approximation for our purpose. 
Here, we use the value of the parameters determined 
by the Ntt reactions in the vacuum, ttia = 1232MeV, 
r = llOMeV, and cta = SOfm^ The medium effects 
on the cross section will be discussed later. Substitut- 
ing TTiN = 940MeV and TOtt = 140MeV, one obtains the 
T dependence of the mean time ta presented in Fig. [T] 
The figure shows that the mean time is ta = 3 ~ 4fm for 
T = 150 ^ 170MeV. One can confirm that the mean 
time hardly changes even for moving nucleons in the 
range of momentum p < 3T by extending Eq. ()14p to 
cases with nonzero nucleon momentum. The lifetime of 
A resonances is rr = l/F ~ 1.8fm. 

The mean time evaluated above is much shorter than 
the lifetime of the hadronic stage in relativistic heavy ion 
collisions. According to a dynamical model analysis for 
collisions at RHIC, nucleons in the hadron phase con- 
tinue to interact for a couple of tens of fm on average at 
midrapidity ^2§\. As a result, at the RHIC energy each 
nucleon in a fireball has chances to undergo the charge 
exchange reactions several times in the hadronic stage. 

Two remarks are in order here. First, the above result 
on the time scales shows that the reactions to produce A 
proceed even after chemical freezeout. These reactions do 
not contradict the success of the statistical model, which 
describes the chemical freezeout [1^, because chemical 
freezeout is a concept to describe ratios of particle abun- 
dances such as (Np) I (Np) and the above reactions do 
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not alter the average abundances in the final state. The 
success of the model, on the other hand, indicates that 
creations and annihilations of (anti-)nucleons hardly oc- 
cur after chemical freezeout. Second, we note that the 
dynamical model in Ref. [l^ uses an equation of states 
having a first order phase transition in the hydrodynamic 
simulations for the time evolution above the critical tem- 
perature Tc- Recently, dynamical simulations have been 
carried out with more realistic equations of states ob- 
tained by lattice QCD simulations js^]- The lifetime of 
hadronic stage evaluated in these studies is more relevant 
to this argument. We, however, note that the qualitative 
behavior of the time evolution seems not sensitive to the 
difference in equations of states [30l |. 

While Ntt reactions frequently take place even below 
the chemical freezeout temperature, Tchom, NN annihi- 
latios and productions almost terminate at Tchom- This 
is necessary for the success of the thermal model. For 
Ec.m. — T the cross section of the NN pair annihila- 
tion is largest among all NN and NN reactions. If nu- 
cleons and anti-nucleons are distributed without correla- 
tion, therefore, all NN and NN reactions cease to take 
place at Tchom- This conclusion is, of course, obtained 
also by evaluating the mean time for each reaction us- 
ing the cross sections [131 as in Eq. ((Ti| . After chem- 
ical freezeout, the only inelastic reactions nucleons go 
through are thus Eqs. (fTO| and pT|) . and after each reac- 
tion the nucleon loses its initial isospin information. Only 
after repeating the reactions Eq. (|12l) twice, the ratio be- 
comes 41 : 40, which is almost even. If medium effects on 
the formations and decays of A are negligible, therefore, 
irrespective of the nucleon distribution at the chemical 
freezeout, the isospin of nucleons at the kinetic freeze- 
out can be regarded random and uncorrelated. On the 
other hand, the nucleon number distribution can have a 
deviation from the Boltzmann distribution reflecting the 
dynamical history of fireballs. 

Because of the absence of correlations between isospins 
of nucleons in the final state, once TVn {N-^) nucleons 
(anti-nucleons) exist in a phase space in the final state, 
their isospin distribution is simply given by the binomial 
one. This conclusion leads to the factorization Eq. ([2]) for 
proton and neutron number distribution in the final state 
for an arbitrary phase space. In particular, the final state 
proton and anti-proton number distribution is written as 

g{Np,Np)= J2 V^{Np,Nn,Np,N^) 

= J'{N^,N^)Br{Np;N^)Br{Np;N^). 

(16) 

Unlike in the simple example in Sec. Ill A[ the nucleon dis- 
tribution function T{N-^, Nfq) in this case is determined 
by the time evolution of fireballs and is not necessary of a 
thermal or separable form as in Eq. ([S]); no specific form 
for T{N-^, Nyi) is assumed here or will be assumed in the 
analyses in Sec. IIIII What we have used here is the fact 



that the time scale for the exchange of isospins between 
nucleons and pions is sufficiently short compared to the 
lifetime of hadronic stage after the chemical freezeout. 
On the other hand, the time scale for the variation of a 
conserved charge in a phase space depends on the form 
of the phase space, and can become arbitrary long by 
increasing the spatial volume. When the time scale is 
long, the information of the physics of the early stages is 
encoded in J-{Nn, N^)- 



3. Medium effects 

Next, let us inspect the possibility of medium effects 
on the formation and decay rates of A. In medium, the 
decay rate of A acquires the statistical factor, 

{I- f{E^)){l + n{E^)), (17) 

where f{E) = {e^^-^''^^/^ + is the Fermi distribution 
function and E^ and E^^ are the energies of the nucleon 
and pion produced by the decay, respectively. The first 
term in Eq. I|17p represents the Pauli blocking effect. At 
the RHIC energy, since the Boltzmann approximation 
is well applied to nucleons, the Pauli blocking effect is 
suppressed. The Bose factor (1 -t-n(i?7r)) in Eq. P?|) . on 
the other hand, has a non-negligible contribution since 
TO^ ~ Tchom- As long as all u^Ett) for the three isospin 
states of the pion are the same, however, this factor does 
not alter the branching ratios Eqs. (fT^ and while 
the factor enhances the decay of A. A possible origin for 
the variation of n(-E,r) is the isospin density of nucleon 
number; since the isospin density is locally conserved, 
the isospin density of pions is affected by the nucleon 
isospin. This effect on n(i?7r) is, however, well suppressed 
since the density of pions is much larger than that of 
nucleons below Tchom- Another possible source which 
gives rise to a different pion distribution is the event-by- 
event fluctuation of the isospin density in the phase space 
at the hadronization. It is, however, expected that the 
effect is well suppressed, again because of the large pion 
density. One, therefore, can conclude that the medium 
effect hardly changes the branching ratios Eqs. and 
P^ . The same conclusion also applies to the formation 
rate of A, since the medium effect on the probabilities of 
a nucleon to undergo reactions Eqs. ([9]) - (fTTj) depends 
only on n{E-^). After all, all medium effects on the ratios 
Eqs. dni) and are negligible. 

When the system has a nonzero isospin density, prob- 
abilities Eqs. ([T^ and receive modifications because 
the three isospin states of the pion are not equally dis- 
tributed, although this effect is not large as will be shown 
in Sec. IIII Dl Even in this case, however, the only mod- 
ification to the above conclusion is to replace the prob- 
abilities r and r with appropriate values, since the reac- 
tions Eqs. (O - pTjl still act to randomize the nucleon 
isospins with the modified probabilities determined by 
the detailed balance condition. 



6 



Here, we emphasize that the large pion density in 
the hadronic medium is responsible for the validity of 
Eq. (fT6|) in the final state. In the hadronic medium, there 
are so many pions which can be regarded as a heat bath 
when the nucleon sector is concerned, while nucleons are 
so dilutely distributed that they do not feel other ones' 
existence. 

So far, we have limited our attention to reactions me- 
diated by A(1232). Interactions of nucleons with other 
mesons, however, can also take place in the hadronic 
medium, while they are much less dominant. It is also 
possible that A interacts with thermal pions to form an- 
other resonance before its decay [3l| . All these reactions 
with thermal particles, however, proceed with certain 
probabilities determined by the isospin SU(2) symmetry 
as long as they are caused by the strong interaction. Each 
reaction of a nucleon thus makes its isospin random, and 
act to realize the factorization Eq. 

4. Low beam- energy region 

The factorization Eq. (ITBl) is fully established for the 
RHIC energy. At very low beam energy, however, pions 
are not produced enough and the duration of the hadron 
phase below Tchcm becomes shorter. Nucleons, therefore, 
will not undergo sufficient charge exchange reactions be- 
low Tchem- When the reactions hardly occur, the isospin 
correlations generated at the hadronization remain until 
the final state. At low beam energy, also the density of 
the nucleon becomes comparable to that of pions, and 
pions can no longer be regarded as a heat bath to ab- 
sorb isospin fluctuations of nucleons. The requirements 
to justify the factorization Eq. (fT6)) . therefore, eventually 
breaks down as the beam energy is decreased. This would 
happen when Tchom ^ ^-n, since the abundance of pions 
is responsible for all of the above conditions. From the 
y'SNN dependence of the chemical freezeout line on the 
T-/XB plane [IBl , the factorization Eq. (fT6|) should be well- 
satisfied in the range of beam energy y^snn ^ lOGeV. 

C. Strange baryons 

So far, we have limited our attention to nucleons. Since 
baryons in the final state in heavy ion collisions are dom- 
inated by nucleons, the nucleon number, which is not 
a conserved charge, is qualitatively identified with the 
baryon one. It is, however, important to recognize the 
difference between these two fluctuation observables es- 
pecially in considering higher-order cumulants. The dif- 
ference predominantly comes from strange baryons A and 
S. In this subsection, we argue a practical method to 
include the effect of these degrees of freedom in our fac- 
torization formula. 

Strange baryons produced in the hadronic medium de- 
cay via the weak or electromagnetic interaction outside 
the fireball. A decays via the weak interaction into p7r~ 



and nTr" with the branching ratio 

Pk-^p : Pk-^u ^ 16 : 9. (18) 
On the other hand, branching ratio of I]+ is 

i's+^p : P^+^n ^ 13 : 12, (19) 

while E~ always decays into mi^ . YP decays into A 
via the electromagnetic interaction and then decays with 
Eq. (dU [131. If the A and E multiplets are created with 
an equal probability, the production ratio of p and n from 
their decays is given by 

Pa,e^p : -Pa^s-.™ ^ 9 : 11. (20) 

Actually, because of the mass splitting between A and 
the E triplets, 5m ~ Tchom/2, the production of the E 
triplets are a bit suppressed compared to that of A. This 
makes the above ratio even closer to even. If one can as- 
sume that the correlations between strange baryons emit- 
ted from the fireball are negligible, therefore, the number 
of nucleons produced by these decays can be incorpo- 
rated into Np and 7V„ in Eq. ^ . The nucleon number in 
Eq. then, is promoted to that of baryons. The same 
argument holds also for A and E. 

In short, by simply counting all protons observed by 
detectors in the event-by-event analysis, A'n and N-^ in 
Eq. ([2) are automatically promoted to the baryon and 
anti-baryon numbers, respectively. 

III. RELATING BARYON AND PROTON 
NUMBER CUMULANTS 

In this section, we focus on the cumulants of the baryon 
and proton numbers, and derive formulas to relate these 
cumulants on the basis of the factorization Eq. ([2]) . With 
these relations the cumulants of the baryon number, 
which is a conserved charge, are calculated from experi- 
mentally observed proton number ones. 

In this section, we change the variables in the proba- 
bility distribution function in Eq. ([2]) as 

V{Np, Np- Nb,N^) = VNiNp, Nn, Np, N^), (21) 

where we have replaced the neutron numbers with the 
baryon ones, A^b = Np + Nn and N^ — Np + Nn. It is 
understood that the prescription discussed in Sec. Ill CI is 
adopted for A, E, and their antiparticles. 

A. Probability distribution functions 

Before deriving formulas to relate the baryon and pro- 
ton number cumulants, in this subsection we first remark 
that the distribution functions of these degrees of freedom 
satisfy a linear relation under the factorization Eq. 
This relation explains why the baryon number cumulants 
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can be represented by the proton number cumulants and 
vice versa. 

Let us start with the final state proton and anti-proton 
number distribution function, Eq. (|f 6p . 



g{Np,Np)= J2 V{Np,Np;NB,N^) 

= J2 HNb,N^)M{Np,Np;Nb,N^) 



JVb.A'b 



(22) 



with 



M{Np, Np; Nb,N^) = B,(iVp; Nb)B,{Np; Nq). (23) 



Equation ((22)) shows that the distribution func- 
tions Q{Np,Np) and J^{Nb,N^) satisfy a linear re- 
lation. Since M{Np, Np; Nb, N-g) has the inverse, 
M-^{Nb,Nb]Np,Np), J^{Nb,Nb) is given in terms of 
giNp,Np) as 



nNB,NB)^ J2 g{Np,Np)M-\NB,NB;Np,Np) 



(24) 



The specific form of M~^{Nb, Nb, Np, Np) is easily ob- 
tained by using the fact that the matrix Eq. (1^51) has a tri- 
angular structure, in the sense that M{Np, Np; Nb, Nb) 
takes nonzero values only for Np < Nb and Np < iVg. 
Using Eq. (l24l). the baryon number distribution function 
T{Nb,Nb) [S^I is in principle determined by g{Np,Np). 
In practice, however, this analysis does not work effi- 
ciently since the elements of M~^{Nb, Nb', Np, Np) are 
rapidly oscillating, which results in large errorbars in 
J-{Nb,Nb) determined in this way. In the following, in- 
stead of the distribution functions themselves, we con- 
centrate on the cumulants of J-{Nb,Nb) and G{Np, Np). 



B. Generating functions and Cumulants 

The moments and cumulants of a distribution func- 
tion are defined in terms of their generating functions. 
The moment generating function for the proton and anti- 
proton numbers with the probability V{Np, Np; Nb, N-q) 
is given by 



G{e,e)= V{Np,Np;Nb,Nb) 

Np,Np,Nb,Nb 



(25) 



and the corresponding cumulant generating function 
reads 



K{e,e) = \ogG{e, 



(26) 



Derivatives of Eq. (I^Sj) give moments of 
V{Np,Np;Nb,Nb), 



{n;np 



Qn gn. 



:G{0,9) 



(27) 



as long as the sum in Eq. (1251) converges, while cumulants 
of the proton and anti-proton numbers are defined with 
Eq. daS]) as 



{iSNp)"{SNp) 



Qn gm 
96*" 89"' 



K{e,o) 



(28) 



The first-order cumulant is the expectation value of the 
operator 



{6Np),^{Np), {SNp),^{Np), 



(29) 



while the second- and third-order cumulants are moments 
of fluctuations, such as. 



{SNpSNp)c = (SNpSNp), 



(30) 



and so forth, with SNx ^ Nx - (Nx)- 

Substituting the explicit form of V{Np, Np; Nb,Nb) in 
Eq. ((ini) for K{e,e), one obtains 

K{e,9)^log J'iNB,NB)eMkNB,N^{9,e)), (31) 
where 



k 



log Y Br{Np; NB)e''-' + log ^ Br{Np; iVg)e 



Ns 



(32) 



is the cumulant generating function for two independent 
binomial distribution functions. With Eq. ([32]). one eas- 
ily finds that this function satisfies fc7VB,-/Ve(0, 0) = and 



^A:jVe,iV3(0,0) = Cn^B, 
gm 

-^JVb,7V3(0,0) -Cm^B, 
^JVb,JV3(0,0) = 0, 



Qn-\-m 



(33) 
(34) 
(35) 



=61=0 



for positive integers n and m, with the cumulants of the 
binomial distribution function normalized by the total 
number 

ii=r, a^r{l-r), 6 - r(l - r)(l - 2r), 

^4 =r(l-r)(l-6r + 6r2), (36) 

and the same formulas for the anti-particle sector. Im- 
posing Eqs. (pij) - (1551) as the structure of K{6,6), cu- 
mulants of net proton and baryon numbers, iVp'^"*'' = 

Np — Np and A''g"°*'' = A^b — Nb , respectively, are calcu- 
lated to be 
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and 



+ SmSNB + ^2SN^)^) + ma^NB - ^^SN^K^Nb ~ 6^^b)) + {^iNb + ^Nb), 



^/ 3el~66 3g-Ci6 \ 



+ 12 ((^'er^jiVp - ^i^^'sNp) i^^'sNp - er'^iVp)) + 3 ((6cr'<57Vp + ^2^^'sNpy 



\ a 



I5$l-i0666 + ef6 
^7 



(37) 
(38) 
(39) 

(40) 



(41) 
(42) 

(43) 



(44) 



A detailed description of the procedure to obtain these 
results is given in Appendix [Xj We emphasize that no 
exphcit form of T{Nb, -ZVg) is assumed in deriving these 
resuhs. Moreover, in Appendixl^we only use Eq. (pij) for 
the structure of K{d, 9) and Eqs. ([55)) - ([55]) for properties 
of kNB,N-s {S, 0) to derive Eqs. (|37)) - p4|) . Therefore, these 
results hold for any distribution functions satisfying these 
conditions with the appropriate choice for the values of 
6 and 6- 

C. Isospin symmetric case 

In hot medium produced by heavy ion collisions, (anti- 
)proton and (anti-)neutron number densities are in gen- 



eral different because of the isospin asymmetry of col- 
liding heavy nuclei. In relativistic heavy ion collisions 
at sufficiently large ^snn and small impact parameters, 
however, the isospin density is negligibly small because a 
large number of particles having nonzero isospin charges 
(mainly pions) are created and most of the initial isospin 
density is absorbed by these degrees of freedom (see, Ap- 
pendix[B]). When the isospin density vanishes, r and f are 
to be set at 1 /2 in the binomial distribution functions in 
Eq. Substituting 

6 = ^, 6 = J, 6=0, 6 = -^, (45) 

into Eqs. ([57]) - (pij) . which are obtained by putting r — 
1/2 in Eq. (l36l) . one obtains 



(4 



not)\ 



(net; 



5N, 



(tot)> 



16 



(46) 
(47) 
(48) 
(49) 
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and 



)^2^ 



=4((<5iV( 
=8((<5iV( 



net)-j2\j 
nct)^3^ 



2(Ar^*°'^), 

I 



12(((5Af(*°'))2) -26(Af(*° 



(50) 
(51) 
(52) 
(53) 



which are the results given in Ref. [25|. Here a note 
is in order about the terms on RHSs of Eqs. ([5T|) - ([5^ . 
Each term on RHS of these equations is not necessar- 
ily uncorrelated with each other. In particular, generally 
J^{Nb,Nj^) is not separable, i.e., it cannot be written as 
J^(A^b,-^b) = fi^B)g{N^). If there is such correlation, 
the statistical fluctuations of these terms are not inde- 
pendent but mutually correlated. Thus, an appropriate 
care need to be taken in estimating the statistical error 
for LHSs of Eqs. (t5T|)-(p|). 



where C and D are constants determined by the chemical 
freezeout condition such as the volume of the system, the 
rapidity coverage, and so on. These relations lead to 



(Nn) (Np) 



(55) 



and thereby r — I — r. One thus can parametrize r and 
f as 



D. Effect of nonzero isospin density 

As the collision energy is lowered, the effect of nonzero 
isospin density eventually gives rise to non- negligible con- 
tribution to the above relations. To investigate this ef- 
fect, we first assume that the isospins of nucleons, anti- 
nucleons, and pions in the final state are in chemical equi- 
librium, as is indicated by the fast Ntt reactions discussed 
in the previous Section. Because the nucleon distribution 
is well approximated by the Boltzmann distribution, the 
numbers of (anti-)protons and (anti-)neutrons in the fi- 
nal state are given with the isospin chemical potential /ii 
and the temperature T as 

{Np) = Ce«/(2^) , {Np} = i?e-w/(2^) , 

(7V„) = Ce-w/(2T)^ ^^^^ ^ i?eW/(2^), (54) 



1 _ 1 

r = --a, r = - +a, 



(56) 



with the negative isospin density per nucleon 



1 {N„) - {Np) _ 1 I - ew/^ 

2 ■ {Nn) + {Np) 2 ■ 1 + ew/T' ■ 



(57) 



a assumes a positive value in heavy ion collisions. 

When the value of a is small, a ^ 1, the effects of 
nonzero isospin density on Eqs. (pij) - are well de- 
scribed by the Taylor series with respect to a. Substitut- 
ing Eq. (|56p in these equations, up to the first order in a 
Eqs. dMl) - (ESI) become 



(nct)\ 



{(SN^^'^f) =4(((5Ar^"^*))2) + 2(iV^*°')) +4a {4{SN^''^'^SN^'°'^) - 3{N^' 

+ 4a (l2((JArj"^*))25iV^t°t)) - 18((JA^j"^*))2) - 6{{SNj,'°''^)^) + 13{N^'°'>, 
((5A^^"'*V)c =16((5A^^"^*))4)c - 48(((5iV^"°*))2,5Ar^'°*)) 



(58) 
(59) 

(60) 



+ 48((<5A^^"^*))2) + 12((,5iV^*°'))2) - 26(iV^*°*)) 
'4a (32((5A^^"'=*))3,5iV^t°t)), - 72{{SN^'"''^f) - 48(5A^^"'=*)(5A^^'°*))^) + l64{6N^'^''''^SNj,' 
-75(7V^"°*M +0(a2). 



(61) 



Next, let us estimate the value of a in relativistic heavy ion collisions. Under the chemical equilibrium condition. 
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the ratio of the charged pion numbers, {N^+) and {N^^- 
having isospin charges ±1, is given by 



-2mi/T 



(62) 



where we have adopted Boltzmann statistics for pions, 
since the effect of Bose-Einstein correlation on the pion 
density is about 10% for Tchcm = "t-tt and does not af- 
fect our quahtative conclusion. The experimental result 
for (A^7r-)/(A^7r+) in the final state is almost unity for 
high energy collisions in accordance with the approxi- 
mate isospin symmetry. Substituting Eq. (|62p in Eq. (157 
and using {N^-) / {Nj^+) — 1 <C 1, one obtains 



1 f {N^-) 
""8 \{N^^) 



(63) 



The value of a, as well as {N.,^-) / {N^+) — 1, grows as 
y^SNN is lowered. In order to see how these parame- 
ters become non-negligible for small -^snn, we focus on 
the 40GeV collision at the SPS ~ 9GeV). For 

this collision, the experimental value of (N^-) / {N^+) is 
1.05 ± 0.05 [29]. Substituting the worst value within la, 
{N„-)/{N^+) = 1.1, in Eq. one obtains a ~ 1/80. 

On the other hand, below the top SPS energy the pro- 
duction of anti-nucleons is well suppressed and one can 



replace all (57V^"°*^ and SN^^""' in Eqs. dH 



(tot) 



6Np to a good approximation. Equation ([61 
pie, then becomes 



dHj) with 
for exam- 



(('5<"V)c 



16(l + 8a)(((5iVp)4) 



48{l + Wa){{dNpf) 
26(1 + 11.5a){Np). 



60(1 + 10.1a)(((5iVp 



(64) 



This result shows that for a = 1/80 the corrections of 
nonzero isospin density to Eqs. ([5(1)) - (15^ are less than 
10% in magnitude. The effect is smaller in relations for 
the lower-order cumulants, Eqs. ([58| - (l60l) . and formulas 
for proton number cumulants, Eqs. ([^5]) - ([^5)1 . 

With these results, one can conclude that the formu- 
las for the isospin symmetric case, Eqs. (|50[) - ()53p . can 
safely be used to the analysis of the baryon number cu- 
mulants for y'SNN — 9GeV with a precision of less than 
10%. Because the production of isospin charged particles 
increases as ^snn goes up, the value of a, and hence the 
effect of nonzero isospin density on Eqs. ([SO)) - are 
more suppressed for higher energy collisions. 

As y^snn is lowered, the value of a grows and even- 
tually approaches the one in the colliding heavy nuclei, 
aA — 0.1. For a ~ 0.1, the first-order correction in 
Eq. ([64]) is comparable with the zeroth-order one. Re- 
lations for the isospin symmetric case, Eqs. (1501) - (I53p . 
therefore, are no longer applicable. For such collision en- 
ergies, however, conditions required for the factorization 
Eq. ([2]) themselves break down as discussed in Sec. Ill Bl 

Before closing this subsection, we recapitulate that the 
suppression of the isospin density in the nucleon sector. 



and hence a, in the final state is caused by the production 
of the large number of particles having isospin charges, 
especially charged pions. In Appendix |BJ we present an 
analysis for this effect. 



IV. DISCUSSIONS 



Recent experimental results on proton number 
cumulants 



As emphasized in the previous sections, the cumulants 
of the proton and baryon numbers are in general differ- 
ent. One, therefore, has to be careful when comparing 
theoretical predictions on baryon number cumulants with 
experimental proton number ones. In this subsection, we 
show that the deviation from the thermal distribution in 
baryon number cumulants becomes difficult to measure 
in proton number cumulants using relations obtained in 
the previous section with some additional assumptions. 

In general, it is possible that, while the net baryon 
number fluctuations in the final state have a considerable 
deviation from the grand canonical ones reflecting the 
hysteresis of fireballs and/or the global charge conserva- 
tion, baryon and anti-baryon numbers separately follow 
the thermal (Boltzmann) distributions. For example, if 
the net baryon number fluctuations above Tc survive un- 
til the final state, the net baryon number fluctuations re- 
main small compared to the thermal ones in the hadronic 
medium, while baryon and anti-baryon number fluctua- 
tions separately follow the thermal one. Generally, cu- 
mulants of net numbers cannot take arbitrary values; for 
instance, the second-order cumulant is constrained by the 
Cauchy-Schwartz inequality: 



< ((<5<"'='^)2) 



(65) 



The values of net baryon number cumulants satisfying 
these constraints are not forbidden. Suppose that, as an 
extreme case, the net baryon number fluctuations com- 
pletely vanish and the left equality in Eq. ()65l) is realized. 
A baryon and anti-baryon distribution function 



(66) 



which is a constrained baryon and anti-baryon number 
distribution following the canonical distribution, con- 
stitutes such an example. The distribution function 
J-{Nb,Nj^) for free gas in the grand canonical ensemble, 
i.e., an unconstrained case, on the other hand, is given 
by Eq. ([8]). 

Now, let us consider the difference between the net 
baryon and net proton number cumulants when the 
baryon and anti-baryon number distributions follow 
Boltzmann statistics while the net baryon number does 
not. Because of the Boltzmann nature of A^b and Nq, dis- 
tributions of Np and Np are also poissonian from Eq. ^ . 
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Thus, cumulants of the baryon and proton numbers sat- 
isfy 

= 2((57Vp)2)HG - 2((<5iVp)3)HG - 2((,5iVp)4),3G, (67) 

and the same for the anti-baryon and anti-proton num- 
bers, where (■)hg is the expectation value for free hadron 
gas (HG) composed of mesons and nucleons at Tchem, i-e., 
a simplified version of the HRG model [s^. The factors 
two in front of the proton number cumulants in Eq. ((67)) 
are understood from Eq. (|6]). 

Using Eq. the second terms in Eqs. (gT]) and pg)) 
are transformed as 

(iv(*°')) = 2{i6N,f + (<57V,)2)hg 

= 2((57V("^*))2)hg, (68) 

= 2{{SNpf - {SNpfUo 

= 2((<5iV(-*))3)HG, (69) 

where in the last equalities we have used the fact that the 
proton and anti-proton numbers do not have correlations 
in the free gas, i.e., {6NpSNp)uG = {iSNp)^6Np)uG = 
{SNp{SNp)^)HG = 0. Substituting Eqs. §EI and §^ in 
Eqs. (|47| and p8)) . respectively, one obtains 

((5Ar(nct))2^ ^ i((<5<°*V) + i((5iV("°'))')HG, (70) 

((57V(-*))3) = i((<5<°*V) + J((<5iV(-'))3)HG. (71) 

These results show that the second terms on the RHSs, 
which come from the binomial distributions of the nu- 
cleon isospin, have significant contribution to the cu- 
mulants of the proton number, and the contribution of 
the net baryon number cumulants, (((^A^g"'^''')"), are rel- 
atively suppressed. Since the second terms give the ther- 
mal fluctuations, these results show that the deviation of 
{{SN^'^*^^)") from the thermal value is hard to be seen 
in the proton number cumulants. Although one cannot 
transform the fourth-order relation Eq. (j49l) to a simple 
form as in Eqs. (170)) and (|7T)) . from the factor 1/16 in 

front of {{6n'^''^^)-^)c in Eq. ^ it is obvious that the 
direct contribution of this term to experimentally mea- 
sured ((<5iV^'"'^)4), is more suppressed compared to the 
lower-order cumulants, and that its experimental con- 
firmation is more difficult. These analyses strongly in- 
dicate that, even if the baryon number cumulants have 
considerable deviation from the thermal values, they are 
obscured in the experimentally measured proton num- 
ber cumulants due to the redistribution in isospin space. 
Such a tendency seems to become more prominent for 
higher-order cumulants. It is known that higher-order 
cumulants of the baryon number have large critical ex- 
ponents and thus can have significant enhancement in 
the vicinity of the critical point [l^. The above result. 



however, indicates that such enhancement is suppressed 
by a factor 1/2" and difficult to measure in experiments 
in proton number cumulants. The analysis of the baryon 
number cumulants with Eqs. ()50l) - ([5^ enables to re- 
move the thermal contribution in the proton number cu- 
mulants and makes the direct experimental observation 
of signals in (((5A'^"°''')")c possible. 

The ^snn dependences of proton number cumulants 
are recently measured by STAR collaboration at RHIC 
[1, [1] . The experimental result shows that ratios between 
net proton number cumulants follow the prediction of the 
HRG model within about 10% precision. We, however, 
emphasize that one should not conclude from this result 
that baryon number cumulants also follow the prediction 
of the HRG model within 10% precision. As demon- 
strated above, the binomial nature of isospin distribu- 
tion makes proton number cumulants close to the ones 
in the HRG model. In this sense, it is interesting that the 
experimental results for skewness and kurtosis neverthe- 
less have small but significant deviations from the HRG 
predictions @. The deviation, for example, in skewness, 

can be a consequence of ((5iV^"°'V) in Eq. dni), which 
possibly reflects the properties of the matter in the early 
stage. 

A remark on Eqs. (|7D)) and ()7ip is in order. These for- 
mulas are obtained with the assumption that baryon and 
anti-baryon number distributions are poissonian, while 
the net baryon number is not. When one further as- 
sumes that the net baryon number cumulants also follow 
the thermal distribution in these results, these formulas 
simply reproduce the free gas result 

{(SNt'^r), = 2(((5A^^"^*))"),,HG (72) 

as they should do. This is easily checked by substituting 

((,5iVi"'='y) = (((5ivi"°*^)")HG in Eqs. m and dZIl)- 

In this subsection, we considered the experimental re- 
sults on proton number cumulants using the results in 
Sec. mil More direct application of these formulas, i.e., 
to determine baryon number cumulants from experimen- 
tal results on proton number cumulants with Eqs. ()50p 
- ()53p . is to be done. The baryon number cumulants 
obtained in this way are to be compared with various 
theoretical predictions. 

B. Efficiency and acceptance corrections 

So far, we have considered the reconstruction of the 
missing information for the neutron number in experi- 
ments using Eq. ([2]). It is possible to extend this argu- 
ment to infer different information on the event- by- event 
analysis. 

An example is the evaluation of the effect of efficiency 
and acceptance of detectors. The experimental detectors 
usually do not have 27r acceptance. Moreover, protons 
entering a detector are identified with some efficiency. If 
one can assume that protons (anti-protons) in the final 
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state is detected by the detector with a fixed probabil- 
ity a (tj) independent of momentum, multiplicity, and so 
on, and the efficiency for each particle does not have cor- 
relations, the distribution function (7V^°'''\ iV^"*""^) 
for the observed proton and anti-proton numbers, Np°^'''^ 

and Np°^^\ respectively, are related to the one for all 
particles entering the detector, Np and Np, as 



(obs) N 



g{Np, Np)B,{N^°'°^^;Np)B^{N^°'"'^;Np), (73) 



or substituting this result in Eq. ([2]) and using the prop- 
erty of the binomial distribution one obtains 

p(obs)(^^obs)^ ^^obs). 

(74) 



E 



Eq. ([71)) indicates that when the deviations of a and a 
from the unity become large, they affect cumulants with 
different orders differently. The effect of efficiency, there- 
fore, cannot be canceled out by taking the ratio between 
cumulants. In particular, as a and a become smaller, 
g(obs)^^^obs)^ ^^obs)^ approach the product of indepen- 
dent Poisson distributions irrespective of the form of 
J^{N^, Nf^). This would be another reason of the present 
experimental results on proton number cumulants 
which is consistent with the HRG model. 

Other experimental artifacts which have not taken into 
account yet in experimental analyses are background 
and misidentified protons. In particular, according to 
Ref. j34i], the contamination from knockout protons is 
not negligible. By their nature, they give poissonian con- 
tribution and make observed proton number cumulants 
approach the poissonian values. Indeed, the HIJING -I- 
GEANT simulation in Ref. ^ shows that these effects 
are considerable. 



of the nucleon and anti-nucleon number distributions in 
isospin space, which is justified for y/suN ^ lOGeV. Al- 
though these results are obtained for isospin symmetric 
medium, the effect of nonzero isospin density in relativis- 
tic heavy ion collisions is well suppressed in this energy 
range because of the abundance of the created pions. 
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is supported in part by Grants-in-Aid for Scientific Re- 
search by Monbu-Kagakusyo of Japan (No. 21740182 and 
23540307). 



Appendix A: Baryon and proton number cumulants 

In this Appendix, we derive Eqs. (|37)) - (|44p . To obtain 
these relations, we start from the cumulant generating 
function Eq. dSU, 



if(0,0) = log^exp [kie,t 



(Al) 



where '^p is a shorthand notation for 
Satb Afg ^b)- this Appendix, we also 

suppress the subscript in k]\j^^N^{6,d). 

We require the following four conditions for the prop- 
erties of k{d, 6): 



k{0,0) =0, 
^fc(0,0) =6.A^B, 
£lfc(0,0) ^Cn^B, 
-fc(0,0) =0, 



89" 86' 



(A2) 
(A3) 

(A4) 
(A5) 



V. SUMMARY 

The most important results of the present paper is 
summarized in Eqs. (ge]) - (|49]) and Eqs. - which 
are formulas relating baryon and proton number cumu- 
lants in the final state in heavy ion collisions. The baryon 
number cumulants are a conserved charge, and one of the 
fluctuation observables which is most widely analyzed by 
theoretical studies. Our results enable to determine the 
baryon number cumulants with experimental results in 
heavy ion collisions, and hence make the direct compari- 
son between theoretical predictions and experiments pos- 
sible. Such a comparison will provide significant informa- 
tion on the QCD phase diagram. The results Eqs. pS)) 
- I|53p are obtained on the basis of the binomial nature 



for positive integers n and m. Eq. (|A2[) is satisfied for 
probability distribution functions normalized to unity. 
Eqs. dMl) - ^EB are Eqs. ([221) - (133 in the text. Ah 
calculations in this Appendix are based only on these 
constraints on K{9, 9). 



1. Net proton number cumulants 

Using K(9,9), the net proton number cumulants are 
given by 

((5iV("'='))")c= (^^-^j KiO,0). (A6) 
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To proceed the calculation of Eq. (jA6p , it is convenient 
to use the cumulant expansion of Eq. (jAip 

711 F 

=i + Y,k{o,o) + \Y.^5k{e,e)f 



(A7) 



Each term on the far right hand side defines each cumu- 
lant, '^p[k[9, ^)]™, up to the fourth order, with 

5k{e,e)^k{e,e)-^k{e,e), (as) 

F 

J2mo, m - Y.^5k{6, eyf - 3 (j^iskie, 9))A . 

F F \ F J 

(A9) 

Because of Eq. (|A2[) . all k{6, 6) and 6k{9, 9) in a term in 
Eq. (|A7|) must receive at least one differentiation so that 
the term gives nonzero contribution to Eq. (|A6[) . This 
immediately means that the m-th order term in Eq. (|A7p 
can affect Eq. (jA6p only if m < n. 

The first-order net-proton number cumulant, Eq. p7p . 
is calculated to be 

(AT^-t)) ^(a, _ a,-)/^(o,o) - ^(a, - 9,-)fc(o,o) 



F 



{^iNb - 6^b) = {^iNb - Ci^b), (AlO) 



with de = d/d9 and 5g = 8/ 89. In the third equality 
in Eq. (|A10|) . we have used Eqs. (|A3|) and (jAi]) . The 
second-order relation, Eq. is obtained as follows: 

= ^(a, - 9,)2fc(0, 0) + 1 - 8ef i5k{0, 0)f 



= Yidl + dl)k{Q, 0) + 2 X - V [(a, - 9,)<5fc(0, 0)] 



F 



2 ^ 



(All) 



To obtain the third line, we have used Eqs. (lASp and (jA2p 
for the first and second terms, respectively. The factor 
two in the second term comes from the number of the 
outcomes of the application of the two derivatives to the 
two 5k{9,9) in the second line. Eqs. (|A3p and (|A4p are 
used in the last equality. 

Similar manipulations lead to Eqs. and (1^0]) . We 
note that the relation, 



{do - 8,f Y.^5k{9, 9))t = 4! ^ [{8e - 8e)k{9, 9)] \ , 

-F F 

(A12) 



F 



makes the calculation for the fourth-order cumulant more 
concise. 

2. Net baryon number cumulants 

To obtain Eqs. (|4T|) - (j44|) . we start from the following 
relation for the net baryon number cumulants. 



(A13) 



with 8(^ = (,^^8e — ^^^8q. We suppress arguments in 
if(0, 0) and A;(0, 0) throughout this subsection. 

The manipulation of Eq. (lAlSp for 7i = 1 is trivial. For 
n = 2, Eq. (|A13p is calculated to be 

F ^ F 

= - E - 5(1)^ - 9(2)^ 



?1 ?1 ?1 ?1 



In the second line, we introduced a symbol, 
and used the relation. 



(A14) 



(A15) 



(A16) 



where we have used Eqs. (|A3p - (lASp . The last equality 
in Eq. (IA14p comes from the definition of K. 

To proceed to ti > 3, we first introduce the following 
notation. 



8 



Xn,m) = T^9e + (-l)"+'"+^(* ^ 5), (A17) 



a, 



(Ti,m,/) 



T^e + (-1) 



n+m+i+2 



(* ^ *), (A18) 



for positive integers n, to, and ^. 5^ ... „^jfor i > 3 
is also defined as in Eqs. (jAlSp . (|A17p . and (|A18I) . One 
easily finds i) (?(„ ,„ ... /) are invariant under the permu- 
tations of the subscripts, for example, 8(n,m.i) — d(^rn.n,i)i 
and ii) when a subscript is one, it can be eliminated, e.g., 
d{n.m,i) = t^(ri,m)j whilc (9(i) — 8^. With this notation, 
derivatives oi 5k are written as 



8^5k = 8(n)Sk, 

8(n)8(rn)Sk = 8(^n,m,2)Sk, 
8{n)8(m)8{l)6k = 8(^n^rn,l,3)6k, 



(A19) 
(A20) 
(A21) 
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and so forth. 

Using these relations, for example, Eq. (IA13P for ?i = 3 
is calculated as 



F 

= OfK - 3 J2i9pk){d^Sk) - dik 

F F 

= dl)K - 3^(9(2)5A;)(9(i)<5fc) - ^ 9(3)^: 

F F 

= df,)K - 3 (9(2) 9(1) X - 9(2,2) i^) - 9(3) if, 

(A22) 

which leads to Eq. P5|) . In the second and last equalities, 
we used 

diK = ^ {d^Sk f + 3 Y,idpk){d^Sk) + J2 9lk, 

F F F 

(A23) 

d(n)d{m)K = d(^„^„i,2)K + ^(9(„)5A:)(9(„,)(5/c). (A24) 



statistics, since this approximation does not alter the 
qualitative conclusion in this Appendix. Under these 
assumptions, the ratios between the numbers of (anti- 
)protons and (anti-)neutrons in a phase space are given 
in terms of /ij and T as 



Nn Np 



l + 2a 



(Bl) 



with a = Np/{Np + Nn), and the ratio of the numbers of 
tt"^ and 7r~ is given by 



ILL = e^^i/"^. 



(B2) 



With these relations, the total isospin in the phase space 
is calculated to be 



Ni =i(iVp - iV„ - Np + Nn) + N^+ - iV,- 



1 -4a2 



(B3) 



A similar manipulation for n = A leads to 

((<57V("'=")4), ^d^.^K - 69(2)9(\)X + 129(2,2)9(i)X 
+ 3df^)K - 49(3)9(i)if - 159(2,2.2)if 

(A25) 



109(2,3)if-9(4)if, 



which gives Eq. 



Appendix B: Isospin density in final state 



with the number of charged pions -/Virch — -^7r+ + N^^- . 

In the initial state of heavy ion collisions, the isospin 
asymmetry of the colliding heavy nuclei aA is approxi- 
mately (Nn - Np)/{2{Np + Nn)) ~ 0.1. Assuming that 
this isospin asymmetry equally distributes along the ra- 
pidity direction in the final state, one has Ni/N^'"''^ ~ 
a A- With Eq. ()B3|) . one then obtains 




aA- 



(B4) 



In this Appendix, we demonstrate that the isospin den- 
sity of nucleons in the final state of heavy ion collisions is 
suppressed owing to the abundant production of particles 
having nonzero isospin charges. 

To simplify the calculation, we consider a gas com- 
posed of nucleons and pions in chemical equilibrium, and 
assume that pions and (anti-)nucleons obey Boltzmann 



The term in the parentheses is larger than unity, and 
becomes larger as more charged pions and anti-nucleons 
are produced. Equation (jB4p thus shows that the value 
of a is more suppressed than aA owing to the production 
of these particles. If the contribution of other particles 
with nonzero isospin charges is taken into account, the 
value of a is further suppressed. 
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